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We obtain the phase diagram and thermodynamic behavior of the Kondo necklace model for 
arbitrary dimensions d using a representation for the localized and conduction electrons in terms of 
local Kondo singlet and triplet operators. A decoupling scheme on the double time Green's functions 
yields the dispersion relation for the excitations of the system. We show that in d > 3 there is an 
antiferromagnetically ordered state at finite temperatures terminating at a quantum critical point 
(QCP). In 2-d, long range magnetic order occurs only at T = 0. The line of Neel transitions for 
d > 2 varies with the distance to the quantum critical point QCP \g\ as, Tn oc \g\^ where the shift 
exponent ip = l/(d — 1). In the paramagnetic side of the phase diagram, the spin gap behaves as 
A ~ \f\g~\ f° r d — 3 consistent with the value z = 1 found for the dynamical critical exponent. We 
also find in this region a power law temperature dependence in the specific heat for faT S> A and 
along the non-Fermi liquid trajectory. For ksT <C A, in the so-called Kondo spin liquid phase, the 
thermodynamic behavior is dominated by an exponential temperature dependence. 



I. INTRODUCTION 

Quantum phase transitions from an antiferromagnetic 
(AF) ordered state to a non-magnetic Fermi liquid in 
heavy fermions (HF) systems have been the subject of 
intense research recently^*. In contrast to classical phase 
transitions, driven by temperature, quantum phase tran- 
sitions can be driven by magnetic field, external pres- 
sure or doping. The physics of heavy fermion com- 
pounds is mainly due to the competition of two main ef- 
fects: the Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teraction between the magnetic ions which favors long 
range magnetic order and the Kondo effect which tends to 
screen the moments and produce a non-magnetic ground 
state. These effects are contained in the Kondo Lat- 
tice Model (KLM) Hamiltonian 2 which can be derived 
from the more fundamental Anderson lattice model in 
the case of well-developed local spin moments^. Although 
this model neglects charge fluctuations, since the relevant 
QCP in heavy fermion materials is associated with a mag- 
netic transition, to consider only spin fluctuations^ turns 
out an excellent approximation. The KLM has been 
studied by different techniques and the general physical 
picture that arises is that at T = 0, there is a quantum 
phase transition from a magnetic phase at small coupling 
strength J/t to a non-magnetic dense Kondo phase at a 
critical value of J/t. 

Here we investigate a simplified version of this model, 
the so called Kondo necklace model 2 (KNM), which for 
all purposes can be considered as yielding results similar 
to the original model. We use the bond-operator ap- 
proach introduced by Sachdev et.al£ that was employed 
previously to both KLM 6 and KNM 4 models but always 
at zero temperature. We extend this approach to finite 



temperatures^ and obtain the phase diagram of the KNM 
for arbitrary dimensions. We find, that this method 
yields a critical Neel line in three dimensions but not in 
2-d as expected from general arguments 8 . This implies 
that fluctuations are taken into account to an important 
extent in spite of the mean-field type of decoupling used 
to deal with the bond-operator Hamiltonian. We also 
calculate the thermodynamic properties along the non- 
Fermi liquid (NFL) trajectory above the QCP 9 ' 10 ' 11 ' 12 . 

The KNM replaces the hopping term of the conduction 
electrons by an X — Y interaction among the conduction 
electron spins 2 . It is given by, 



H = t 



E 

<i,3> 



( ^x ^x 



rfr^ + J^.n 



(1) 



where and Si are independent sets of spin-1/2 Pauli 
operators, representing the conduction electron spin and 
localized spin operators, respectively. The sum de- 
notes summation over the nearest-neighbor sites. The 
first term mimics electron propagation and in one dimen- 
sion can be mapped by the Jordan- Wigner transforma- 
tion onto a band of spinless fermions. The second term 
is the magnetic interaction between conduction electrons 
and localized spins Si via the coupling J. 

The paper is organized as follows: In Sec. II we intro- 
duce the bond operator treatment of the KNM. In Sec. 
Ill we make and discuss the approximations required to 
solve the problem. The Green's function method is em- 
ployed to attain this solution, as it can be easily gener- 
alized to finite temperatures. We obtain a closed set of 
equations for the Green's functions that can be imme- 
diately solved without no further approximations. The 
temperature dependence of the singlet and triplet order 
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parameters is found in Sec. IV. In Sec. V the finite tem- 
perature phase diagram of the KNM is discussed. In Sec. 

VI we study the paramagnetic phase and find the depen- 
dence of the spin gap with the distance to QCP. In Sec. 

VII the behavior of the specific heat of the KNM is pre- 
sented. The last section (Sec. VIII) offers the conclusions 
and discussions. 



II. BOND OPERATORS REPRESENTATION 

We use the method of bond operators, which has the 
advantage of making the connection to the lattice degrees 
of freedom most direct. For two S = h spins, Sachdev et. 
al. 5 introduced four creation operators to represent the 
four states in Hilbert space. This basis can be created 
out of the vacuum by singlet \s) and triplet \t a ) = 4,|0) 
(a = x,y, z) operators. In terms of these triplet and 
singlet operators the localized and conduction electrons 
spin operators are given by, 



T i,a = 2 ( _ S i^i,a ~ A,a S i ~ ^ e ap 1 t\ [} t in ). (2) 

where a, (3 and 7 take the values x, y, z, repeated in- 
dices are summed over, and e is the totally antisymmetric 
Levi-Civita tensor. In particular, an important feature of 
the bond operator mean-field approach is that the sim- 
plest mean-field theory to be used below already yields 
ground states and excitations with the correct quantum 
numbers; so a strong fluctuation analysis is not needed 
to capture the proper physics. The restriction that the 
physical states are either singlets or triplets leads to the 
constraint s's + J2 a A^a — 1- Moreover, the singlet and 
triplet operators at each site satisfy bosonic commuta- 
tion relations [s,s^\ — 1, = 5a,pi [ s iAi] — 0- Sub- 
stituting the operator representation of spins defined in 
Eq. (2) into the original Hamiltonian and considering 
the commutation relations, we obtain 

H = Hq + Hi + H2 

where 

i a 

+ Z)w(*i«i+D*U*i,a- 1 ). 

i a 

H l = \ E [ S i s Xa*j,a + 4 s Aa*L + ' 

(3) 



with a = x, y. Hq represents the interaction between 
spins S and t in the site i and the constraint s^s + 
S a ta^a = 1 is implemented through the local chemical 
potentials Hi and H2 are terms associated with the 
hopping. There is still another term in the Hamiltonian 
that consists of three triplets and one singlet operator. 
This is not taken into account since it vanishes in the 
approximation that we use below 4 . 

III. ANTIFERROMAGNETIC PHASE 

The Hamiltonian above, at half filling, i.e., with one 
conduction electron per site, can be simplified using a 
mean-field decoupling of the quartic terms while still 
retaining the relevant physics. The resulting effective 
Hamiltonian H m f with only quadratic operators is suf- 
ficient to describe exactly the quantum phase transition 
from the ordered AF state to the disordered Kondo spin 
liquid at least for d > 3, as discussed below. Besides 
it also yields sensible results in 2-d where no line of fi- 
nite temperature transitions is found. Relying on the 
nature of the strong coupling limit (t/J) — * we take 
(s|) = (sj) = s, which corresponds to a condensation of 
the local Kondo spin singlets on each site . Next to de- 
scribe the condensation of one local Kondo spin triplet 
tk, x on the AF reciprocal vector Q = (jr/a, it /a, n/a), we 
introduce: t-^ x = VNiS^Q + Tjk^ corresponding to fixing 
the orientation of the localized spins along the x direction. 
The quantity t is the mean value of the x-component spin 
triplet in the ground state and 77k, x represents the fluctu- 
ations. Finally the translation invariance of the problem 
implies that we may assume the local chemical potential 
as a global one. 

We will consider here only the terms Hq and Hi; the 
Hi term has only small contributions to the results there- 
fore it is neglecte d 13 i 14 . 

After performing a Fourier transformation of the boson 
operators, the mean-field effective Hamiltonian is given 
by, 

H mf = N^Js 2 +f is 2 -^ 

+ (7 + /*) E 4a,» + (1 + » - 

^ ' k 

+ E [ A ^i,A, x + A k ('/k , k., + »M-k,x) 

k 

+E [ A ^Lv, + A k (AjU, y + «k,»t-k,»)] , 

k 

(4) 

with A k = uj + 2A k , A(k) = £) s =i cosk s , A k = 
jts 2 X(]i), N is the number of lattice sites, Z is the to- 
tal number of the nearest neighbors on the hyper-cubic 
lattice. The wave-vectors k are taken in the first Bril- 
louin zone and the lattice spacing was assumed to be 
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unity. This mean-field Hamiltonian can be solved using 
the Green's functions to obtain the thermal averages of 
the singlet and triplet correlation functions. These are 
given by, 



< t 



k,a ' k,a 



1 



2ir uj 2 — oj^ ' 



t 1 1 

Z7T UJ — UJq 



(•5) 



(0) 



temperatures we calculate in this section the finite tem- 
perature order parameters s 2 = s 2 (T) and t 2 = t 2 (T). 
The free energy can be directly obtained from the en- 
ergy of the excitations given by the poles of the Green's 
functions found in the previous section. It is given by, 



2 N 
G = £o~pJ2 Ml + n(«k)] - - ln[l + n(wo)] (9) 



/3 



where 



where, a = x,y. The poles of the Green's functions de- 
termine the excitation energies of the system as ojq — 
( 4 + /i), which is the dispersionless spectrum of the lon- 



gitudinal spin triplet states and uJk — ±-y/AjJ 



*k - (2A k ) 2 

that correspond to the excitation spectrum of the trans- 
verse spin triplet states for both branches oj x = oj y . From 
these modes and their bosonic character an expression for 
the average energy at finite temperatures can be easily 
obtained, 



SO 



^ Wk ( coth - 1 



(7) 



where 



£o = N 



/is 2 



1 - 2 

- + M - -tZs 2 



(8) 



is the ground state energy of the system and (3 = X/ksT, 
s and t the singlet and triplet order parameter respec- 
tively. Since the parameter s is always nonzero^ and 
t 7^ in the antifcrromagnetic phase, we minimize the 
ground state energy with respect to t to find, /i = ^Zts 2 — 
J/4 and consequently, w k = \ZU 2 yJ\ + 2X(k)/Z. The 
ground state energy £o corresponds to a magnetically 
long-range ordered state characterized by the momen- 
tum q = Q. The excitations over this ground state are 
given by two spin-wave branches associated with t x and 
t y that represent rotations or transverse oscillations of 
the order parameter i. The third mode associated with 
t z corresponds to fluctuations in the amplitude of this 
order parameter. In the next section we consider low 
temperatures with the purpose of studying the behavior 
of the system near the magnetic instability. 



IV. SINGLET AND TRIPLET ORDER 
PARAMETERS AT FINITE TEMPERATURES 

The low temperature thermodynamic and transport 
properties of heavy fermions in the vicinity of the mag- 
netic quantum critical point are far from being com- 
pletely understood^. In order to study the KNM at finite 



n H = 2 I coth — - 1 



(10) 



Minimizing the free energy of the KNM using 
(de/dn,de/ds) = (0,0), we can easily get the following 
saddle-point equations, 



1 + 


J 


l 






~zt ~ 


2N 




J 


1 


1 - 






Yt 


2N 



2A(k) coth ^ 
Z 2 



1 



: COth 



(3uj k 



2A(k) 

z 



(11) 



where £ = ^(coth^j 2 - — 1). Generally the equations 
for s and t in Eq. (fTTj) should be solved and at 
T = the results of Ref. (4) are recovered. For 
J/t > (J/<) c , triplet excitations remain gapped and at 
J/t < (J/t) c , the ground state has both condensation of 
singlets and triplets at the antiferromagnetic wave vec- 
tor Q = (tt, it, 7r). Then the quantum critical point at 
T = 0, (J /t) c separates an antiferromagnetic long range 
ordered phase from a gapped spin liquid phase. For fi- 
nite temperatures, the condensation of singlets s occurs 
at a temperature scale which, to a first approximation, 
tracks the exchange J. The energy scale below which the 
triplet excitations condense is given by the Neel temper- 
ature (T/v) which is calculated in the next section. 



V. CRITICAL NEEL LINE IN THE KNM 

In terms of the KNM, the condensation of triplets (sin- 
glets) follows from the RKKY interaction (Kondo effect). 
Thus, the fact that at the mean-field level, both s and t 
do not vanish may be interpreted as the coexistence of 
Kondo screening and antiferromagnetism in the ordered 
phase 4,6 for all values of the ratio J/t < (J/t) c . Notice 
that the bond-operator mean-field theory is appropriate 
in the strong coupling limit and near the QCP. However 
it does not give an accurate description of the weak limit 
for (J/t) — * 0, where the ground state is macroscopically 
degenerate. 

The Neel line giving the finite temperature instabil- 
ity of the antiferromagnetic phase for J/t < (J/t) c is 
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obtained as the line in the T vs (J ft) plane at which t 
vanishes (t = 0). From Eq. (fTT|) we can then obtain the 
boundary of the AF state. From this equation we get, 



theorem^. For d > 3, f(d) is finite and the equation for 
the critical line is given by, 



T N oc \g\ 



(17) 



151 
Z 



2N ^ 

k 



1 



2A(k) 



coth 



-i +e, (12) 



where g - 
the QCP. 

IN Sk 



\{J/t) c — (J/t)\ measures the distance to 
The latter is given by, (J/t) c = Z[l — 
where Z is the number of nearest 



^l+2A(k)/Z 

neighbors. Expanding close to the wave- vector Q = 
(ir, 7r, 7r) associated with the antiferromagnetic instabil- 
ity we get, 



A(k) 



y + 0(fc 4 ), 



(13) 



This yields the spectrum of transverse spin triplet exci- 
tations as, 



U>o\ 1 + 



2A(k) 
Z 



Dk. 



(14) 



where yD = Wo/v 2d, with loq the z-polarized dispersion- 
less branch of excitations and d is the Euclidean dimen- 
sion. Notice that the low temperature specific heat in 
the antiferromagnetic phase is dominated by the contri- 
bution of these modes. In 3d this has a T 3 temperature 
dependence due to the linear, phonon-like, dispersion re- 
lation. 

Replacing Eq. (fl4|) in Eq. (fT2|) and considering that 
for temperatures fc^T <C loq, £ goes to zero faster than 
the first term of Eq. (fT2|) we obtain, 



l£l 
Z 



2 d - 2 S d Z d / 2 fk B T 



LJ 



d-1 



,d-2 



(cothw — 1) du, (15) 



where Sd is the solid angle and u = f3uiQk/2\/Z. For 
temperatures k B T <C loq the integral in Eq. (| 1 5(1 can be 
calculated and we get 



jd-2 



s d z d " 



Z TT^" 1 



(16) 



where f(d) = J °° u d ~ 2 (coth u — 1) du. The solution of 
this expression gives us the critical line of Neel transi- 
tions for any dimension and temperatures k B T <C ujq 
where ujq = {i + /i) is the dispersionless spectrum of 
the longitudinal spin triplet states that tracks J. We 
notice that the integral f(d) diverges for d < 3 showing 
that there is no critical line in two dimensions at finite 
temperatures^, in agreement with the Mermin- Wagner 



which defines a shift exponent ij> = l/(d — 1). In par- 
ticular for 3d, f(3) — 7r 2 /12 and the critical line of Neel 
is given by k B T^ = O.lSwolffl 1 ^ 2 - In summary, in this 
section we have obtained analytically the expression for 
the Neel line close to the QCP and we have shown that 
this line does not exist for d = 2 as expected. 



VI. PARAMAGNETIC PHASE 

Our interest here is to obtain the thermodynamic prop- 
erties near, but above the QCP, i.e., for finite tempera- 
tures in the spin liquid phase for J/t > (J/i) c . In this 
case we have to consider only the condensation of singlets 
since t = in the paramagnetic phase. We find, 



— r I coth 



(3ujq 
2 



coth 



(18) 



The ground state energy is now given by e' Q — 



N(-ps 2 + 
ergy by, 



M s - W + 



Sk ( w k — wq) and the free en- 



G' 



N 

n(wk)] - — ln[l 



■n(w )] (19) 



The transverse spin triplet excitation spectrum can be ex- 
pressed in the general form Wk = \/A£ — (2Ak) 2 , where 
the excitations are triply degenerate magnons (triplons) 
with a spin gap that vanishes on approaching the QCP, 
as we show below. 

We minimize the paramagnetic KNM free energy, de- 
riving the following saddle-point equations: 



2-s 2 



N Z-~i v ; 9 



(3u> k 



N 



Wk 



1 E- 



Ak , , /3w k 
, / — coth — — 
N ^— ' w k 2 

k 



2f • (20) 



At T — the self-consistent equations above are equiv- 
alent to those found in the Kondo spin liquid ground 
stated. Considering again the expansion close to the QCP 
of Eq. (|13p . we obtain the transverse spin triplet excita- 
tion spectrum in the paramagnetic phase as, 



UJy = y/A 2 + Dk 2 



(21) 



These quasi-particle excitations form a band whose band- 
width is a function of t. The band-minimum is at k = tt 



5 



and the triplet excitation spectrum has a spin gap given 
by, 



wo\ 1 



yz 

2 ' 



(22) 



In Eq. [2TJ L> = yio 2 /2 and y = is 2 /wo is a dimension- 
less parameter. The spin gap decreases as the hopping t 
(hence y) increases and vanishes when y — > 2/Z. 



A. Dependence of the spin gap on the distance to 
the QCP 

We now obtain the dependence of the spin gap A with 
the distance to QCP, given by g = \ J/t — (J/t) c \, in the 
quantum disordered Kondo spin liquid state from Eq. 
([20]) at T = 0. This was done previously for the spin gap 
in spin ladder systems but always in one dimensio n 13 ' 16 . 

The first step is to express Eq. ([20]) as an equation for 
the parameter y, 



1 



2 1 - d K d VI + yd 



K(rj) 



with 



tt/2 



sin 



(23) 



(24) 



where r\ = \/2yd/ (yd + 1) and K (if) in one dimension 
is the complete elliptic integral of the first kind^ 3 -. In 
Eq. (|2"2"|) we see that when y = y c = 2/Z, the spin gap 
vanishes yielding 



2 

Vc 



1 



TT d VI + Vcd 



K(Vc) 



(25) 



where rj c = v / 2y c d/(y c d + 1). Since \g\ = \ J/t - (J/t) c 
we obtain, 



where C(Vc, d) = ^^TTTP^WI^ with C(^,2) « 
297.01 and C(j/c, 3) ~ 17.78 in two and three dimensions 
respectively. /i/J and y are given by the self-consistent 
equations (|2"P|) and (l2"3"j) . respectively and are never crit- 
ical. This relation between the spin gap and \g\ shows 
that when (J/t) decreases from its strong coupling limit, 
the triplet gap at the wave vector k — (ir, it, it) decreases, 
and vanishes at \g\ = 0. Since A cx Iffl 1 / 2 , when \g\ — > 0, 
we can immediately identify the gap exponent vz = 1/2 
at the quantum critical point of the Kondo lattice. 



VII. SPECIFIC HEAT IN THE 
PARAMAGNETIC PHASE 



The purpose of this section is to consider the param- 
agnetic half-filled KNM energy given by Eq. ([19)) and 
discuss its finite-temperature properties near the QCP. 
Already the ground state and finite-temperature prop- 
erties of the KLM have been intensively studied for 
one dimension using the density matrix renormalization 
group (DMRG) metho d 17 ' 18 ' 19 . For two dimensions this 
has been accomplished using finite-temperature Lanczos 
technique^, perturbation theory^ and quantum Monte 
Carlo (QMC) simulation^. Lately, also for the KNM us- 
ing QMC for two dimensions- 3 . However, less analytical 
work has been carried out for higher-dimensional KNM 
and KLM. This way, motivated by the above, we are 
going to analyze the specific heat C /T = 7 using an ana- 
lytical treatment valid for any dimension. We obtain 7 in 
the paramagnetic disordered phase, in two cases: First in 
the so called non-Fermi liquid trajectory (\g\ = 0, T 7^ 0) 
and then, below the spin gap for fc^T <C A cx Iffl 1 ^ 2 - We 
consider the free energy Eq. (fT9"]) with ui^ = VA 2 + Dk 2 
and take into account only the contribution of the trans- 
verse spin triplet excitation spectrum cjk, since the con- 
tribution of longitudinal triplet excitation tends to 
zero faster at low temperatures ksT < wo ~ J. 



where 
G(y) 



2(y c - y) . 

y c y TT a 



1 



yVl + yd 



K(v) 



1 



y c Vl + y c d 



K(Vc), (27) 



The above equation can be expanded in Taylor's series 
for y near y c yielding, 

G(y) « |- 2—^ K(r,)\ y = yB (y-yc)+0(y-y c ) 2 , (28) 
oy z/Vi + yd 

Now considering Eq. (j26[) we can express the spin gap 
energy by, 



A = J 



yz\g\ 



4 J) V 2(Z + t(y c ,d)Y 



(29) 



A. Specific heat along the NFL trajectory 

Non Fermi-liquid behavior is often found near a mag- 
netic QCP—, indicating that the NFL state in those sys- 
tems may be linked to the magnetic instability at T = 0. 
In this state the coefficient of the linear term of the spe- 
cific heat does not saturate as expected from the Landau 
scenario but shows a temperature dependence as the tem- 
perature is lowered. Exactly at the QCP, this may occur 
down to the lowest temperatures 2 ^. Next, we are going 
to find analytically C/T at the NFL trajectory. From 
Eq. (0 and using C/T = -d 2 G'/dT 2 we get, 



C/T 



D d / 2 2n 



- d s(d)(k B T) d -\ 



(30) 



6 



where 



u 3 (u 2 



4) 



' (coth 2 — — l)du, 



with x q = A/k B T , A = Lj Q ^/yZ\g\/[2(Z + ((y c ,d))], 
u = yjxl + D/3 2 k 2 and g = \{J/T) - (J/T) c \ measures 
the distance to the QCP. Equation (|30|) yields the expres- 
sion for the specific heat in the paramagnetic region and 
it is due, as already we pointed out, to the contribution of 
bosons t x and t y . At \g\ =0 and low temperatures, such 
that, J 3> k B T ^> A (NFL trajectory) we can rewrite 
Eq. ((201) as follows, 



C / T \\g\=o - %|=o 



S d k 2 B ,(\g\ = ,d) i 

£>d/2 2n d 



{k B Ty-\ (si) 



where q{\g\ = 0, d) = J u d+1 (coth 2 § - l)du is 24C(3) 
and 167r 4 /15 in two and three dimensions, respectively 
and £ is the Riemann zeta-function. Thereby C/T oc 
T d ~ x at the QCP. Notice that this is consistent with 
the general scaling result C/T oc with the dy- 



namic exponent taking the valui 



,25 



1 



Together 



with the previous result for the gap exponent vz = 1/2, 
we find that v — 1/2, a result in agreement with the 
mean-field or Gaussian character of the approximations 
we have used to deal with the bond-operator Hamilto- 
nian. Furthermore, since z = 1, in three dimensions 
d e f f = d+ z = d c = A where d c is the upper critical 
dimension for the magnetic transition^. Consequently, 
the present approach yields the correct description of the 
quantum critical point of the KNM for d > 3. 



B. Specific heat in the Kondo spin liquid state 

In the region II of the Fig. (1) we have to consider the 
paramagnetic contribution to the specific heat given by 
Eq. (fT5|) but now taking into account that k B T <C A. 
Then for k B T < A w y/\g\ we obtain, 



C/T = 



2S d A 



d+2 



n d D d/2 keT 3 



„A/k B T 



6(d,T), 



(32) 



where cosh it/ = 1 + x Q 
/ U/ (sinh w) d_1 cosh 3 udu. 
we obtain S(d = 2,T) = 
3,T) 



V2 
3 



2(^)3/2 + 17 ( 



x = A/k B T and 5(d) = 
This can be calculated and 
;igr + and S(d = 

17 ( k B T \b/2\ 



A 



? ) in two and three 



dimensions respectively. In this case the specific heat is 
governed by the exponential term and we can conclude 
that the dimensionality does not make much difference in 
thermodynamic properties. This point was also indicated 
using numerical methods for the 2c? KLM2£. 
Summarizing we have 



f T d ~ l 
C/T cx < e -A/k B T 



NFL trajectory, J > k B T > A 
, k B T -C A 



(33) 



NFL Trajectory 

T a|ql 1/(d " 1) r d " 1 
n |y L y a T 



PARA 




Aa|g| 



1/2 



•' y a e 



A/k „T 



|g|=o 



J/t 



FIG. 1: (Color online) Schematic phase diagram of the Kondo 
necklace model at finite temperatures. The AF phase is lo- 
cated below the Neel line Tjv (solid line), which vanishes at 
a critical value (J/t) c (QCP). The spin gap energy (dashed 
line) is shown as a function of g. Below this line there is a 
Kondo spin liquid state. The temperature dependence of the 
specific heat 7 = C/T is given along the NFL trajectory and 
in the Kondo spin liquid state for any dimension. The phase 
diagram shows that at low temperatures a small value of J/t 
supports the development of AF order via the RKKY interac- 
tion. If the value of J/t becomes too high, the local moments 
are quenched by the conduction electrons and the system has 
a dense Kondo ground state. 



where in the latter case only the dominant term has been 
written. Then, we have obtained the specific heat due 
to the magnetic degrees of freedom in the paramagnetic 
phase of the half-filled KNM for any dimension and tem- 
peratures k B T <g; J. We have investigated it in two cases: 
along the NFL trajectory and in the dense Kondo spin 
liquid phase where an exponential temperature due to the 
spin gap dominates over the power law associated with 
different dimensionalities. Unfortunately, there arc few 
existing results in higher dimension considering the ap- 
proaches done; almost all results above were obtained for 
the strong coupling limit (t = 0). Our work should pro- 
vide a valuable benchmark for other approaches. While 
we have not attempted a detailed comparison on or fit to 
experiment, this would be possible. 



VIII. CONCLUSIONS 

In the present work, we have constructed analytical ex- 
pressions to find the phase diagram of the Kondo neck- 
lace model for any dimension and at low temperatures, 
such that, k B T <§; J. Although there were several ap- 
proaches treating the Kondo lattice model with methods 
similar to the one we use, they were all restricted to zero 
temperature. The present finite temperature treatment 
allowed us to determine the critical line for antiferro- 
magnetic transitions, as well as, the thermodynamic be- 
havior along the quantum critical trajectory and in the 
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non-critical part of the phase diagram. We obtained the 
relevant critical exponents governing the critical line and 
the thermodynamic behavior due to the magnetic degrees 
of freedom. For three dimensions these results turn out 
to be exact since d e f f = d + z = A is equal to the upper 
critical dimension for the Kondo lattice. 

We also found that the spin gap in the Kondo singlet 
phase vanishes as \/\g\ at the QCP consistent with z = 1 
and a mean field correlation length v = 1/2. 
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